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Abstract 

Backlund-Darboux transformations are closely related to the in- 
tegrability and symmetry problems. For the generalized Backlund- 
Darboux transformation (GBDT), we consider conservation laws, ra¬ 
tional extensions and bispectrality. We use the case of the nonlinear 
optics equation (and its auxiliary linear system) as an example. 


1 Introduction 

Backlund-Darboux transformations and related commutation methods con¬ 
stitute one of the most fruitful approaches to the construction of explicit 
solutions of linear and nonlinear equations and to the explicit spectral theo¬ 
retic results (see, e.g., [HI^IHISHTn] and numerous references therein). Similar 
to various symmetries studied in group analysis, Backlund-Darboux transfor¬ 
mations (BDTs) transform solutions of linear differential equations of some 
hxed class into solutions of differential equations of the same (or another 
hxed) class and potentials into potentials, whereas potentials are often solu¬ 
tions of the corresponding integrable nonlinear equations. Some connections 
between Darboux transformations and potential symmetries are discussed in 
an interesting paper na. 

Though the idea of BDTs and commutation methods goes back to Back- 
lund, Darboux and Jacobi, the notion of the so called Darboux matrix has 
a much later origin. According to J. Cieslihski [Tj, ’’all approaches to the 
construction of Darboux matrices originate in the dressing method”. In 
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order to explain the notion of the Darbonx matrix, we consider some initial 
linear system 

Ux = G{x,z)u, n G C™', (1.1) 

where ■= -^u {u^ = in the case of several variables), G{x, z) is some 
m X m matrix fnnction (i.e., G{x, z) G C™'^™') and C stands for the complex 
plain. The Darboux matrix (matrix fnnction) w satishes the eqnation 

Wx = Gw-wG, G{x,z) (1.2) 

It is easy to see that if fll.ip and fll.2p hold, the prodnct 


u = wu 


satishes the transformed system 

Ux = G{x, z)u. 

The cases, where the Darbonx matrix w can be constrncted explicitly, if u 
and G are given explicitly, are of special interest. Clearly, in those cases we 
constrnct explicitly u and also G. 

We are interested in a generalized version of the Backlnnd-Darbonx trans¬ 
formation for which we nse the acronym GBDT. In GBDT (for the case of one 
space variable) Darbonx matrix is presented as the transfer matrix fnnction 
in Lev Sakhnovich form: 

(1.3) 

where 7415 — SA 2 = LIiII^, (1.4) 

Im is the m X m identity matrix. We note that an additional variable x 
appears in wa, Bi, 112 and S', so that we deal with matrix fnnctions nA;(x), 
S{x) and Darbonx matrix wa{x, z). GBDT is applicable to all systems Ux = 
G{x, z)u, where G rationally depends on z, that is, 

Gix, 2 :) = ( ^ z^Qkix) + - Cs)~^qskix) 

k=0 s=l k=l 

(see the paper [13], the reviews [T3lll5| and Ghapter 7 in [T6| as well as 
nnmerons references therein). The coefficients % and Qsk of the transformed 
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(rational) matrix function G{x, z) are expressed via the coefficients and 
Qsk of G{x, z) and matrix functions S{x), ni(a:) and 112 (x). Matrix functions 
nfc(x) are constructed as solutions of systems related to the initial system. 
Here, we shall consider in detail GBDT for the system 

Ux = G{x,z)u, G = izD — [D,g\, D = dia.g{di,d 2 , ■ ■ ■ ,dm} = D*; 

( 1 . 5 ) 

B ^ diag{5i,52,... ,5m} ih = ±1); [D, q] Dq - qD, 

( 1 . 6 ) 

where diag denotes a diagonal matrix. System fll.Sp is an auxiliary system 


for the well-known iV-wave (nonlinear optics) equation: 

[D, gt] - [D, Qx] = [[D, p], [D, p]], (1.7) 

D = diag{di,...,dm} = D*. (1.8) 

More precisely, equation fll.7p is equivalent to the compatibility condition 

Gt-Fx + [G,F]=0 (1.9) 

of the auxiliary systems 

Ux = G{x,t, z)u, G = izD — [D,Q], (1-10) 

Ut = F{x,t, z)u, F = izD — [D,g]. (l-H) 


We shall study conservation laws, rational extensions and bispectrality for 
the case of GBDT for system fll.51) and also for the Wwave equation fll.7p . 

The integrable model describing interaction of three wave packages [18] 
is the most well-known subcase of the iV-wave equation fll.7p . Putting, for 
instance, m = 3, B = G and di> d 2 > ds, and using transformations m Gh. 
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di)/(di — d2). 

'4^2 — (da — d2)/(d2 — da). 



^3 = (da - di)/(di - da); 

( 1 . 12 ) 

iy^di — d2 Pi2, 

02 = —W d2 — da P2ai 

03 = —iy/di — da pia, 

( 1 . 13 ) 
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we rewrite fll.7p in the standard form of the corresponding 3-wave interaction: 


(0l)t + (^ 2 )* + '02(02)x — 

{h)t + '03(</>3)x = i£0102, (1-14) 

where 0(a;, t) stands for the function, which takes values complex conjugate 
to (j){x, t), and 

£ =(did2 — <^2^1 “1“ (^2*^3 — <^3(^2 4“ d ^ d \ — did ^') 

X ((di - d2)(di - d3){d2 - (1.15) 

If a is a scalar value or matrix, the notation a stands for the scalar, which 
is complex conjugate to a, or the matrix with entries complex conjugate to 
the entries of a, respectively. 


2 Preliminaries 

GBDT for the system (ll.5h and for the iV-wave equation was described in [13] 
(see also [Ml Subsection 1.1.3 and Section 7.1] and references therein). In 
this section, we give some necessary dehnitions and results from [13] and m 
Subsection 1.1.3 and Section 7.1]. We consider the case a; > 0 for the system 
(II.5p and the case x > 0, f > 0 for the iV-wave equation. We fix n x n 
matrices A and iSlO) = S'(O)* and an n x m matrix 11(0) such that 

71^(0) - ^(0)71* = in(o)5n(o)*. ( 2 . 1 ) 

We introduce the nxm matrix function n(x) and the nxn matrix function 
S'(x) via initial values 11(0) and S'(O) and differential equations 

= -iTinzi + n[zi, p], s^ = udbu*. ( 2 . 2 ) 

Relations (12. Ih and (12.2p yield AS{x) — S{x)A* = in(x)i?n(x)*. Then the 
following proposition describes GBDT determined by 7l, S'(O) and 11(0). 

Proposition 2.1 Let system (II.5p [such that (II.Op holds) be given. Then 
{in the points of invertibility of S{x)) the matrix function 

wa{x, z) = Im- iRn(x)*^(x)-^(7l - zIn)~^A{x) (2.3) 
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is a Darhoux matrix of fll.Sp and satisfies the equation 

z) = {izD - [D, g{x)])wAix, z) - wa{x, z) {izD - [D, ^(x)]), (2.4) 

where 

g= g- BU*S-^n, g* = BgB. (2.5) 

Moreover, if detS'(x) 7 ^ 0 for x > 0, a normalized fundamental solution 
w{x, z) of the transformed system 

Wx = {izD - [D, g\)w, w{0, z) = Im (2.6) 

is given by the equality 

w{x,z) = wa{x, z)w{x, z)wa{0, z)~^, (2.7) 

where w is the fundamental solution of the initial system fll.51) normalized by 
w{0,z) = Im- 

Next, we add the variable t and consider n(a;, t), S{x, t) and wa{x, t, z) which 
are determined by A, 11(0, 0) and 5(0, 0) = 5(0, 0)* via eqnations fl 2 . 2 p . 02.31) 
and 

Hi = -iAUD + n[5, g], St = UDBU*. (2.8) 

Instead of 02.ip we assnme that 

715(0,0) - 5(0,0)71* = in(o, o)sn(o, o)*, ( 2 . 9 ) 

which implies that AS{x,t) — S{x,t)A* = iIl{x,t)BIl{x,t)*. Proposition 12.11 
yields: 

Proposition 2.2 Let an m x m matrix function g {g* = BgB) be contin¬ 
uously differentiable and satisfy the N-wave equation 01.7p . Then g of the 
form 

g{x, t) := g{x, t) — BIi{x, t)*S{x, t)“^n(a;, t) (2-10) 

satisfies {in the points of invertibility of 5) the equality = BfiB and the 
N-wave equation. 

Using m Therem 6.1] on wave fnnctions, we easily obtain the next state¬ 
ment. 
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Remark 2.3 If the conditions of Pro'position \2.^ hold, det S{x,t) ^ d on the 
semi-band 0<x<oo, d <t < e and w{x,t,z) is the initial wave function 
[i.e., Wx = Gw, Wt = Fw and w(0,0,2;) = Im), then the transformed wave 
function w{x,t,z) is given by the equality 

w{x, t, z) = wa{x, t, z)w{x, t, z)wA{d, 0, z)~^. (2-11) 

Example 2.4 The real-valued case is of special interest m Section S.f], 
It is immediate from Proposition \2.^ {and formulas defining g considered in 
this proposition) that equalities 

g = g^ A = -A, n(o,o) = n(o,o), 5(o,o) =:s(o,o), ( 2 . 12 ) 

yield the equality g = g. Assume that (12.121) holds, and sog = g. Then, in the 
subcase m = 3, B = I 3 , di > d 2 > d^ of the 3-wave equation, we can rewrite 
fll.l4p in the real-valued form and obtain solutions of the corresponding exact 
resonance equations. Namely, we set 

Ti = ~i0i = —^/di — d2 gi2, T 2 = ~i02 = —\/d2 — g23, 

T3 = “103 = —Vdi — ds gi3, 

and, taking into account fll.l4p . derive 

{(pi)t + 'ipi{(pi)x = e(p2(p3, {^ 2)1 + i’2{T2)x = 

{T3)t + '4’3{T3)x = —£TiT2 - (2.13) 

3 Conservation laws, rational extensions and 
bispectrality 

3.1 Conservation laws 

When g is differentiable with respect to t (to x), we can differentiate both 
sides of the hrst relation in fl2.2p (in fl2.8p ) with respect to t (to x). If p is 
differentiable with respect to x and t and gt is continuous, then, according 
to the so called Clairaut’s (or Schwarz’s) theorem, we have Axt = B-tx- Thus 
Tlxt = B-tx is the necessary condition of the compatibility of the hrst relations 
in (12.2p and (12.8p . In a similar way, Sxt = Stx is a necessary condition of the 
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compatibility of the second relations in fl2.2p and fl2.8p . In other words, the 
necessary compatibility conditions are 


i-iAUD + U[D, g])t = i-iAUD + U[D, p]),, (3.1) 

(3.2) 

where Ifa, should be substituted by —lAAD + n[Zl, p] and lit should be sub¬ 
stituted by —lAAD + n[Il, p]. 

Proposition 3.1 Let q he differentiable and satisfy the N-wave equation 
cai. Then the equalities fl3.ip and fl3.2l) hold. 

Proof. Substituting expressions for fl^ and fla, we obtain 

{-iAUD + Il[D, p])t = - iA{-iAIlD + n[5, g])D 

[—lAIlD + n[zi, g])[D, p] -|- n[zi, Pi], 
{-iAUD + n[5, g])^ = - iA{-iAUD + U[D, g])D 

(— iyinzi -i- n[zi, g])[D, p] -|- n[zi, p^,]. 

Therefore, taking into account fll.7l) we derive 

{-iAUD + U[D, g])t - {-iAUD + U[D, pj), 

= U{[D, Pi] - [3, p„] - [[D, p], [3, p]]) = 0, 

and so fl3.ll) is proved. Using the same substitutions as before, we have also 

{UDBU*)t = ( - iAu3 + n[5, g])DBU* + UDB{i3u*A* + [p*, 5]n*), 

(3.3) 

{u3bU*)^ = {-iAUD + U[D,g])3BU*+ u3B{iDU*A*+ [gfD]U*). 

(3.4) 

Equality 03.21) follows from 03.31) . 03.41) and Bg* = gB. □ 

The sufficiency of the proved above equalities U^t = Utx and Sxt = Stx 
for the compatibility of 02 .2p and 02 .8p is not self-evident but may be proved 
in a way similar to the proof of [Ml Theorem 6.1]. 

Remark 3.2 ITe note that relations 03.11) and 03.21) may be considered as 
conservation laws for n(a;, t). 
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3.2 Rational extensions 

According to Propositions 12.11 and 12.21 in the case 


p = 0, (j(A) = {0} (a stands for spectrum), (3.5) 

the transformed potentials p are rational matrix functions {rational exten¬ 
sions in the terminology of [H]). More precisely, if p = 0 and A is nilpotent, 
formulas fl2.2p . fl2.3p . fl2.5p . fl2.7p . fl2.8p . fl2.10p and fl2.1ip imply the following 
proposition. 

Proposition 3.3 Let fl3.5p hold. Thenll{x) and S{x) are matrix polynomi¬ 
als with respect to x, n(a;, f) and S{x,t) are matrix polynomials with respect 
to X and t; p(x) = p(x) / det S(x) and g{x, t) = p{x, t) / det S{x,t), where p{x) 
is a matrix polynomial with respect to x, p{x, t) is a matrix polynomial with 
respect to x and t, det S{x) is a polynomial with respect to x and det S{x, t) 
is a polynomial with respect to x and t. Moreover, for Darboux matrices wa 
we have wa{x, z) = P(x, z)/ det S{x) and wa{x, t, z) = P{x, t, z)/ det S{x, t), 
where P{x,z) {P{x,t,z)) is a matrix polynomial with respect to x and 1/z 
(x, t and ^jz'). Thus, the transformed fundamental solutions w{x,z) and 
wave functions w{x, t, z) are expressed via matrix polynomials and exponents 
of diagonal matrices: 

w{x, z) = wa{x, z)e'^^^^WA{^^ 

w{x, t, z) = wa{x, t, z)P^^^^~^^^^wa{0, 0, z)~^. (3.6) 

Proof. Since p = 0, the hrst equation in fl2.2p yields that the ith column fi 
of n is given by fi{x) = e"‘‘^*"^^/i(0) = ^{-idixA)’^fi{0), and so n(x) 

is a matrix polynomial. (Here we used the equality A^ = 0, which holds for 
all n X n nilpotent matrices.) Hence, the second equation in fl2.2p yields that 
S{x) is a matrix polynomial. In the same way we prove the polynomial form 
of n(x,f) and S{x,t), and the required properties of g follow. According 
to formula (1.84) from [16], p. 24] we have wa{z)Bwa{z)* = B and, in 
particular, the equalities 

'ifA(0, z)~^ = Bwa{0, z)*B, tCA(0, 0, z)~^ = Bwa{^-, 0, l)*B 

are valid. Finally, the representations of wa{x, z) and wa{x, t, z) follow from 
the dehnition (12.3p . properties of H and S and the expansion {A — zln)~^ = 


















Remark 3.4 System fll.Sp includes important subclasses: self-adjoint and 
skew-self-adjoint Dirac-type systems {which are also called Zakharov-Shabat 
or AKNS systems). Namely, putting 


and B = 
system 


D=j = 


‘■mi 

0 


-L 


(mi + m 2 = m), V = 


0 V 
V* 0 


(3.7) 


j, we rewrite fll.Sp . fll.6p in the form of the self-adjoint Dirac-type 


Ux = i(^j + jV)u, V = 2i^i2, (3.8) 

where Q 12 is the upper right mi x m 2 block of g. Assuming fl3.7p and B = Im, 
we rewrite fll.Sp . fll.6p in the form of the skew-self-adjoint Dirac-type system 


Ux = {izj + jV)u, V = -2gi2- 


(3.9) 


3.3 Bispectrality 

The notion of bispectrality was introduced in [2] (see also [7] and references 
therein). Matrix bispectrality for system fll.Sp was studied in im. Namely, 
it was shown that for g constructed in Proposition 13.31 and, correspondingly, 
for the solution W{x, z) = wa{x, of system 

W^ = {izD-[D,'^)W, (3.10) 

there is a non-degenerate matrix linear differential operator B with respect 
to the variable 2 ; such that B{z)W = 0. It was suggested in [7] (and in 
some earlier works by F. Griinbaum and coauthors) that the bispectrality 
requirement WB{z) = Q{x)W is more meaningful than the requirement 
B{z)W = Q{x)W. Unfortunately, the approach from [T7] does not work 
properly for that case. Therefore, we should consider the system 

W^ = {izD-[D,'^)W-izWD (3.11) 

instead of system fll.Sp . In view of fl2.4p . the matrix function wa satisfies 
this system in the case p = 0. Thus, we can deal in the same way as in [5], 
where Dirac systems were treated. Namely, if ^{A) is concentrated at one 
point A G C we write down the resolvent {A — zln)~^ in the form 

n 

(A - zh)-^ = - 5^(^ - A)-'=(A - (3.12) 

k=l 

By virtue of fl2.3p and fl3.12p the next proposition is immediate. 
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Proposition 3.5 Assume that g = 0, the spectrum of A is concentrated at 
some point A G C and g is given by fl2.5p . Then system fl3.1ip is bispectral 
in the sense of m- 

Indeed, it is easy to choose some coefficients Cg so that the operator B{z) = 
Cs{z — Xy is non-degenerate and satishes the equalities 

S(^)(^-A)-^ = 0 (l<fc<n). (3.13) 


According to formula fl2.4p and identity p = 0, the matrix function W = 
wa{x,z) satishes fl3.1ip . Therefore, dehnition fl2.3p of wa and equalities 
fl3.12p and fl3.13p imply B{z)W = 0. 

If we want to consider system fl3.10p . where p = 0 and cr(A) = {A}, A G 
M, we rewrite the hrst relation in fl2.2p in the form 

= -i{A - Xlm) (n(x)e‘^^^) D. (3.14) 

Hence, taking into account the fact that (t{A) is concentrated at A, we see that 
n(x)e‘^^'^ = P{x, A) = P{x), where P{x) is a matrix polynomial with respect 
to X. Substituting n(a;) = P{x)e~^^^^ into fl2.2p and using the equalities 
D = D*, X = X, we obtain 

= P{x)e-'^^^^DBy^^^P{xy = P{x)DBP{x)y i.e., ^(x) = P(x), 

(3.15) 

where P{x) a matrix polynomial with respect to x. Introduce an operator 
B = diag{Hi,..., Bm} such that Bk has the form 


N 


= NeN, 


e=i 


and notice that the equality 

'^cekiz)-^ {{z - X)”"wA{x,z)P''^^ek) = 0, 


£=1 


et := (!<<;< m) 


(3.16) 


(3.17) 


is equivalent to the relation 

N ( Q A ^ 

^Qfc(x) l — + idkx] 
e=i ^ ^ ' 

X (^{z - xrp^^^ (j™ - iBP{xyP{x)-\A - zIn)-^P{x)^ = 0, 
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which, in turn, is equivalent to 


N , Q X £ 

i — + idkx\ (3.18) 

i=i ^ ^ ' 

X (^{det P{x)){z - xr [im - iBP{xyP{x)-\A - x4)-'P(x)) e^) = 0. 

Open Problem. Are there some cases, where matrices D, A, S'(O) and 
n(0) generate nontrivial ^ such that WB{z) = 0 for the solution W{x,z) = 
wa{x, of fl3.10p and for P 7 ^ 0 of the form given above? In other 

words, are there D, A, S'(O), n(0) and {cik} satisfying fl3.18p ? 
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